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Abstract 

We study a discrete-time approximation for solutions of systems of decoupled 
forward-backward doubly stochastic differential equations (FBDSDEs). Assuming 
that the coefficients are Lipschitz-continuous, we prove the convergence of the scheme 
when the step of time discretization, |jt| goes to zero. The rate of convergence is ex- 
actly equal to |ti| 1//2 . The proof is based on a generalization of a remarkable result on 
the 2 -regularity of the solution of the backward equation derived by J. Zhang ifTTI . 
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1 Introduction 

In this paper, we study a discrete time approximation scheme for the solution of a system 
of the (decoupled) forward-backward doubly stochastic differential equations (FBDSDEs, 
in short) on the time interval [0, T]: 

X t =x + fob(X s )ds + Jo o(X s )dW s 

(1.1) 

Y t = h(X T ) + Sj f(s,X s ,Y s ,Z s ) ds + f t T g(s,X s ,Y s )dB~ s - f t T Z s dW s . 

Here W and B are two independent Brownian motion such that, the integral with respect 
to B t is a backward Ito integral and the one with respect to W t is a standard forward Ito 
integral. Let us note that such equations naturally appear in probabilistic interpretation of 
stochastic partial differential equations (SPDEs, in short). Indeed, under standard Lipschitz 
assumptions on the coefficients b,G,f,g, and h, the existence and uniqueness of the solution 
(Y,Z) have been proved by Pardoux and Peng [9]. Moreover, they give the link between the 
classical solution of SPDE in the following. More precisely let consider the SPDE 

— — u(t,x) — [Lu(t,x) — f(t,x,u(t,x),a* (x)Vu(t,x))] — g(t,x, u(t,x))()B s = 0, 

u(T,x) = h(x), (1.2) 
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where denotes the Wick product and, thus, indicates that the differential is to understand 
in Ito's sense, and 



i,j 1 J i 1 

Under more strengthen assumptions (the coefficients /, g and h are C 3 class), the compo- 
nent Y of the solution of (jl.ip is related to the classical solution u of SPDE dl.2| ), in the 
sense that 

Y t = u(t,X t ). (1.3) 

Furthermore, Buckdahn and Ma relax the assumptions of coefficient to standard Lipschitz 
one and they proved among other that the relation (jl.3p give the stochastic viscosity solution 
of SPDE (|1.2|) . Thus, solving (|l.ip or (|1.2p is essentially the same. However it is known 
that only a limited number of BDSDE can be solved explicitly. In order to solved the large 
class of BDSDE and of course provide an alternative to classical numerical schemes for a 
large class of SPDE, the numerical method and numerical algorithm is very helpful. 

In the one stochastic case, i.e g = 0, the numerical approximation of (jl.ip has already 
been studied in the literature; see e.g. Zhang [11], Bally and Pages [2], Bouchard and Touzi 
131 or Gobet et al. Q. In Q, the authors suggest the following implicit scheme. Given a 
partition regular grid ji : = to < t\ <....< t n = T of the interval [0, T\, they approximate 
X by its well-know Euler scheme X K and (F,Z), by the discrete-time process (F^,Z")o</< n 
defined backward by 

< 

Y* = E [Y* +1 \? ti ] +A? +1 f(t i ,X*,Y*,Z%), 

where Y* = h(Xj), A n W i+ i = W, i+l - W t , and Af +1 = t i+1 - t { . Then, it turn out that the 
discretization error 

Err n {Y,Z) = \ sup E\Y t -Y t n \ 2 + f E [\Z S -Zf | 2 ] ds 

{0<t<T JO 

is intimately related to the quantity 

£ f " +1 E [\Z S - Z ti \ 2 ] ds where %. = -i-E 

Under Lipschitz continuity conditions on the coefficients, Zhang [ 1 Q was able to prove 
that the latter is of order of |tc[, the partition's mesh. This remarkable result allows them 
to derive the bound Err n (Y,Z) < C|7t| x / 2 . Observe that this rate of convergence cannot be 
improved in general. Consider, for example, the case where X is equal to the Brownian 
motion W, h is the identity, and / = 0. Then, Y = W and Y t f = W u . 
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In this paper, we extend the approach of Bouchard and Touzi Q, and approximate the 
solution of by the following backward scheme. 

z l = a^7 E ? [(C +*fen,*£ 1 ,l£ 1 )A8 m ) AW i+1 ) , 

Y* = Ef [Y* +1 +g(t i+h XZ +1 ,Yl 1 )AB i+1 ] +f(t i ,X*,Y*,Z*)A i+1 

where Y t n = h(Xj) and AB, + i = B tj+l —B tj . By adapting the arguments of Bouchard and 
Touzi [3 ], we first prove that our discretization error Err n (Y,Z) converge to as the step of 
the discretization |7t| tends to 0. We then provide upper bounds on 

max sup E|y f -F,,| 2 + V / E [\Z S - | 2 ] ds. 
<<» 0<t<ti i Ju 

When the coefficients are Lipschitz continuous, we obtain 

n ~ l fti+i 

max sup E\Y,-Y,.\ 2 + V / E \\Z S -Z t .\ 2 ] ds < C\n\. 

i<» 0<t<ti ' i Ju 

This extends to our framework the remarkable result derived by Zhang iTTTTl . It allows us 
to show that our discrete-time scheme achieves, under the standard Lipschitz conditions, a 
rate of convergence exactly equal to Iftl 1 / 2 . 

Observe that, in opposition to algorithms based on the approximation of the Brown- 
ian motion by discrete processes taking a finite number of possible values (see e.g. ifTUll 
and the references therein), our scheme does not provide a fully implementable numerical 
procedure, since it involves the computation of a large number of conditional expectations. 

This paper is organized as follows. In Section 2, we introduce some fundamental knowl- 
edge and assumptions of BDSDEs and give extension of the remarkable L 2 -regularity re- 
sults derived by Zhang [11 J to the doubly stochastic case, which is our first main result. In 
Section 3, we describe the approximation scheme and state convergence result, our second 
main result. 

Notations. We shall denote by M"' the set of all n x d matrices with real coefficients. 
We simply denote R" = M' u and M" = M">". We shall denote by ||a|| = (L-, ; -a 2 j) 1/2 the 
Euclidian norm on M"> d , a* the transpose of a, a k the k-th column of a. To simplify, we 
denote respectively by \x\ and au, the norm and the the fc-th component of a £ R". Finally, 
we denote by x.y = £,x i -;y i - the scalar product in R". 

2 Forward-Backward doubly SDEs 
2.1 Preliminaries and Assumptions 

Let (Hi, jFi,Pi) and (Q.2, ^2,^2) be two complete probability spaces and T > a fixed final 
time. Throughout this paper we consider {W,,0 < t < T} and {B t ,0 <t<T} two mutually 
independent standard Brownian motions processes, with values respectively in M. d and R c , 
defined respectively on (Oi,iFi,Pi) and {Q.2, f 2^2)- For any process (r\ s : < s < T) 
defined on (£2,-, jF ; -,P ( ), (i = 1, 2), we denote 

T?t = o{r\ r -T\ s ,s<r<t}V9{, 7? = 
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In the sequel of the paper unless otherwise specified we denote 

Q. = £li x a 2 , 7 = 7\ <S> Ji and P = Pi ® P 2 . 

Moreover, we put 



7t = 7 t 



where 9£ is the collection of P-null sets and denote F = (7 t )t>o- Further, for random 
variables £((Qi), GOi £ Cli and ^(002), G>2 £ ^2, we view them as random variables in Q. by 
the following identification: 



two functions 



e(cfl) = e(a>i); p(a>) = p(co 2 ), co= ((01,002). 

Given C > 0, we consider two functions b : M. d — > W 1 and a : W 1 — > 
satisfying the Lipschitz condition 

(m)\b(x)-b(x')\ + \\o(x)-o(x')\\ <C\x-x'\, Vx,x' £R d . 

Then it is well-known that (see e.g Karatzas and Shreve (61), for any initial condition x € M. d , 
the forward stochastic differential equation 



Xf = x + f b(X s )ds + f a(X s )dW s , t £ [0, T] 
Jo Jo 



(2.1) 



has a jTf-adapted solution (X t 



t)0<t<T 



satisfying 



E( sup |X f | 
0</<r 



Before introducing the backward doubly SDE, we need to define some additional nota- 
tions. Given some real number p > 2, we denote by S p the set of real valued adapted cadlag 
processes Y such that 



sup \Y t \ p 

0<t<T 



< 00, 



\\Y\\ s p = E 

Lo<, 

tt p is the set of progressively measurable M rf -valued processes Z such that 



E 



\Z t \ p dt 



Vp 



<C 00 ■ 



Mp 



The set <B P = S p x tt p is endowed with the norm 

||(y,z)|| SP = (||F||^ + ||z||^) J 

The aim of this paper is to study a discrete-time approximation of the pair (Y, Z) solution 
on [0, T] of the backward doubly stochastic differential equation 

Z s dW s , 0<t<T. (2.2) 



Y t =h(X T ) + J t f(s,X s ,Y s ,Z s )ds + J g(s,X s ,Y s )dB s - 

By a solution, we mean a triplet (F,Z) G r B p satisfying (|2.2p . 

In order to ensure the existence and uniqueness of a solution to (|2.2[) . and the conver- 
gence of our discrete-time approximation, we assume that the map / : [0, T] x M. d x R x 



t,g: [0,T]x 



and 



satisfied the Lipschitz condition: 
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(H2) 

(0 \f(s,x,y,z)-f( S , ^,y',z')\ z <C{\s- S f+\x-J<f\ 2 + \y-y'\ 2 + \z-e\ 2 ) 
(ii) \g(s,x,y) - g(s,x' ,y')\ 2 <C(\s-s'\ 2 + \x-x>\ 2 +\y-y / \ 2 ) 
(Hi) \h(x)-h(x')\ 2 <C\x-x'\ 2 

for some constant C > independent of all the variables. 

Remark 2.1. In order to ensure the existence and uniqueness to the solution of (|2.2p . we 
need only that / and g are Lipschitz with respect variables y and z. See Pardoux and Peng 
Q for more detail. 

The following lemmas collect without proof, some standard results in SDE and BDSDE 
literature. We list them for ready references. For ease of notation, we shall denote by C p 
a generic constant depending only on p, the constants C, b(0), a(0), h(0) and T and the 
functions /(., 0,0,0) and g(.,0,0). 

Lemma 2.2. Assume b and O satisfy (HI) and X be the unique solution of forward SDE 
(TO) . Then 

\\X\\ P SP < C p {\ + \x\ p ) 

and 

n%-x s \p]<c p {\+\x\p)\t-s\p' 2 . 

Lemma 2.3. Assume (H2) and (Y,Z) be the unique solution of backward doubly SDE (|2.2p . 
Then 

\\(Y,Z)\\l P <C p {\ + \x\n 

and 

E [\Y t - Y S \P\ < C p {(1 + \x\ p )\t ~ sr 1 + ||Z||> ,} . 
2.2 L 2 -regularity 

In this subsection we establish the first main result of this paper, which we shall call the 
L 2 -regularity. Such a regularity, plays a key role for deriving the rate of convergence of our 
numerical scheme in Section 4 and, in our mind generalized Theorem 3.4.3 in iTTTi . 

To begin with, let Ji : = to < ... < t n = T be a partition of the time interval [0, T], with 
|7t| = maxi<;<„ \ti-i —ti\, the size of the partition, and X be the solution of the forward SDE 
(|2.ip . We denote by (F,Z) the solution of the following backward SDE 

Y t = 4>* (X t0 , . . , X tn ) + J f(s,X s , Y S ,Z S ) ds + g(s,X s ,Y s )d%- ^ Z s dW s , (2.3) 

the generalized form of BDSDE (|2.2|) . Next, for X % the well-know Euler scheme of X that 
will be explicit in Section 3, let (Y n ,Z K ) be the adapted solution to the following BDSDE 

Y t K = f{Xl...X) + f f(sXX,Z?)ds+ f g(s,X?,Y?)d% - [ T ZfdW s . (2.4) 
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To simplify presentations, in what follows we assume that X t , Xf G R d , and the other 
processes are all one-dimensional. But the results can be extended to cases with higher- 
dimensional on this processes without significant difficulties. For simplicity we also denote 
by S = (X,Y), © = (X,Y,Z) and E % = {X n ,Y n ), <d % = (X n ,Y n ,Z n ). 

Now we have 

Lemma 2.4. Assume the functions (f : R d (' l+1 ) -> R, / : [0, T] x R d x R 2 -> R and g : 
[0,r]xR d xK-tR satisfying assumptions (H2) wifft adequate norm. For each I <i <n, 
we define 



1 



7^ - - 



/" 



Z^s 



wfere Ef_!(.) = E(.|^ V j/). TTzen 



limsup|7l| E 



max sup \Y t - Y ti l | 2 + £ / |Z, - Z 71 , l 2 ^ 



l <'< n ti^<t<ti 



< oo. (2.5) 



Before prove this important theorem, we state the following needed result. To this 
end let us assume the following: (j)" G C l b {R d ^), f G cJ^flO.T] xl^x R 2 ) and g G 
Cj' 1 ([0, r] x R d x R). Moreover, for all x = (jc , ....,jc„) G R rf (' l+1 ), 



EKWI<C (2.6) 

i=0 

We also design by (p M the partial differential of cp which respect the variable u. 
Next, we denote by VX n the solution of the following variational equation:, 

VXf = I d + f b x (X?)VX?dr + f o x {X?)VX?dW r , (2.7) 
Jo Jo 

and by (VT 71 , V'Z 71 ) the solution of the following BDSDE on [r f _i, T] : 

+ j\g x (Z«)VX? + g y (Z?)V%«)]dl! r -j\ i Z?dW n t G T], 
fon = l,...,n. (2.8) 

On the other hand, we denote by 



v k y k = £V%*l [ti _ uti) (t)+V»Y?-l {T] (t), t G [0,r]; 
i=i 



(2.9) 



hence V 71 }" 1 is a cadlag process. 
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For application convenience, we shall rewrite y K Y K in another form. Note that for each 
i (|2T8|) is linear. Let (7°,^°) and j= 1, ...,n be the adapted solutions of theBDSDEs 

7? = f\fM-W + f y {®«)f+fM)£¥r 

+ ^ [fo(S?)VX? + g, {^ r )f]dB r - ^ t£dW r , (2. 10) 

= ,X?jVX« + j\f y (®*)Y r +f z (Q*)Wdr 

respectively, then we have the following decomposition: 

V% = 7? + £y^ ^[Uft)- (2-11) 

/>* 

We may simplify (|2.1ip further. Let us define, for any r| G L 1 (F, [0, T]) and (®i , ©2) G 

L 2 (F,[0,r];M) xL 2 (F,[0,r];I), 

A?(T]) = exvU\{r)drY s,te[0,T), 
^/(©i) = expj^GiW^-i^l©!^)^^, ^G[0,r], 
2 £/(0 2 ) = exp|^0 2 (ryB r -^'|02W| 2 ^|, S ,fe[0,r]. 

( ! £/(®i ) and 2 E/(@2) are respectively the well known Dalean-Dade stochastic exponential 
of ©1 with respect W and ©2 with respect B). Then it is easily checked that, for any p > 0, 
one has 

Mi&tW = ^©OA'C^^I©,-! 2 ), (2.12) 

and 

f'E/C©,)]- 1 = 'E/C-e.OAfd©.-! 2 ), / = 1,2. (2.i3) 

In particular, we denote, for s,f G [0,T], 

A? = A°(-f y ) 2 ^(- gy ), M* = \'(f z ), (2.14) 

and if there is no danger of confusion, we denote A = A and M = M°. Since f z is uni- 
formly bounded, by Girsanov's Theorem (see, e.g., @) we know that M is a P-martingale 
on [0, T], and W t = W t — J^f z (&^)dr, t G [0, T] is an F-Brownian motion on the new proba- 
bility space (0,5 ,P), where P is defined by S = My. Moreover noting that / y , / z and gy 



7 



are uniformly bounded, by virtue of (|2.12|) and (|2.13|) one can deduce easily from ()2. 14[) 
that, for p > 1 , there exists a constant C p depending only on T, C and p, such that 

E ( sup \At\ p + [A," 1 A < C p ; E ( sup [\M t \P + \M^\ P ]) < C p ; 

\0<t<T J \0<t<T J 

E(\A t -A s \P + \A; l -A; l \P)<C p \t-s\P/ 2 ; (2.15) 

E (\M t -M s \ p + \M~ l -M~ l \ p ) < C p \t - s\ p l 2 . 

Lemma 2.5. Assume o, b E Cl and f, g, I satisfy the previous assumptions. Then for all 
i = 1, ...,n 

y% n = Ut+L^i) M; l A t - fj x {<d*)VX?A; l drA t - £ gx (Z*)VX?A; l d i B r At, 

where ^ and for j = 1 , • • • , n will be explicit in the proof. 
Proof. Let us denote the following: 

1° = fJfA®?)VX?K l dr + J T gx (Z?)VX?A- l d i B n g = ^\ 
^ = f t AT l +flf^)VXfA^dr + ^g^)VXfA-'dB r 

l i = h*(x*,--;X? n )vx? i A- T \ 5 = eV, %=iA7 l - 

Then, using integration by parts and equation (|2.10p we have, for i = 0,1, ....,n, 

%=V-f*'&Wr, te[0,T], 

so that, /q C, l r dW r being a uniformly integrable martingale with in particular zero expectation, 
we get 

? f =E(f'> f ). 

Therefore, by the Bayes rule (see e.g, |'6] Lemma 3.5.3) we have for t € [0, T] 

7? = f t At- f U®*)VX?A7 l drA t - f g^)VX?A- l dB r A t 
Jo Jo 

= $Mi%- f f x {®*)VX?A- l drA t - f ' g x {Z*)VX?A- l d*B r A u 
Jo Jo 

i = %A t = I Af = E (MtI%) A t = QM^At, 

where, for i = 0,1, ....,n, 

%\ = E (MtI%) = E (MtV) + j\[dW s . (2.16) 
8 



Note that the boundedness of/, and (l2~T5j) imply thatM^ G L P {£1) and VX GZ/(F,C([0,r];M rf )) 
for all /? > 2. Therefore for each > 1, ()2.6[) leads to 



E <^ V \M T tJ\ \ < CE\ \M t \ p sup \VX t \ p } < C. 
L=o J I o< f <r 



In particular, for each j = 0, • • -,n, MjQ G L(Tt)- So (|2.16p makes sens. Finally the result 
follows by (IXTTT) . □ 

Proof of Lemma 2.4. For all 1 <i<n and each t G ,?,•), applying Lemma l2"31 we get 

E (|y f -y fi i | 2 ) <c\n\. 

Then by Burkolder-Davis-Gundy inequality we have 



E 



max sup \Y t — Y tj _, 



KKn 



<CJC 



(2.17) 



The estimate for the second term of the left hand in (|2.5p is little involved. First we 
assume that b,a, (j) 71 , /, g G such that tf 1 satisfied (f2T6|) . Let recall (y^Z 71 ) denote the 
adapted solution to the BDSDE (|2.4|) and X 71 the solution of the Euler scheme associated to 
EDS (|2.1|) . Under the Lipschitz conditions on b and a, we have 



lim max E 

7i-*0 KKn 



sup \X?-X t \ 2 + sup l^-X^, 



0. 



(2.18) 



0<t<T ti-i<t<ti 

Now by the Lipschitz assumption on (j) 71 and (|2.18|) . applying Lemma 2.2 we know that 



limEi sup \Y t n -Y t \ 2 + f \Zf - Z t \ 2 dt 1 = 0. 
[o<t<T Jo J 

Recalling (|2.5p and applying Lemma 3.4.2 of Zhang flTTl we have 



(2.19) 



\Z s -Z*J 2 ds 



ti-i 



< 



\Z S -Zf \ 2 ds 

i ((—1 i 



< 



2 £ E [ f (|Z, - Zfl 2 + |Zf - 2* , | 2 )^ 

1=1 L-^fi— 1 



(2.20) 



By ()2.19l) and (|2.20|) . to estimate the second term and prove the theorem it remain to show 
that 



ti-i 



\Z?-Z*J 2 ds 



<cM, 



(2.21) 



where C is independent of 7t. 

To do this, let us recall that from Proposition 2.3 of (3 and its proof, we know that the 
martingale part Z 71 has a continuous version given by 



Zf = VX[VX?]-'o{X?), Vf G [tt-uti) 
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which together with Lemma 2.5 provide 



Therefore, 



izr-z^Kit+it+if+it 



(2.22) 



where 



;>i ;>f,-i + i 



xlM-iVJVX-J^a^J, 



£+E# 



;>«' 



JO 



Recalling ()2. 15|) and applying Lemma l2T2l and Lemma l2~3l one can easily prove that 

lE(|/, 3 | 2 + |/,T)<C|7i|. 
Recalling ()2. 16|) and (|2.6|) . we have 

|S,° + £S/|<CE( sup \VX?\\T t ). 

j>i 0<t<T 

Thus by using again Lemma l2T2l and Lemma l231 one can similarly show that 

E(|/ 2 [ 2 ) <C\n\. 
It remains to estimate l} . To this end we denote 

r,= sup {i + \x?\ + \[vx?]- 1 \ + \m; 1 \}. 

0<s<t 

Noting that A is bounded and that r tj l G F ti l , by ()2. 16|) . we have 



(2.23) 



(2.24) 



E|//| 2 < CE^rJ_, 



fiT+E^-RL+E&U 



< ce jr^E 1 1§? - | 2 + £ 15/ - ^ | 2 |^, 



< ce < r! 



f/-1 







2 




/" 






) 




/>«' 
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Therefore, by following the step of ifTTTl . we get 

£E ^ $\ 2 <ftj < c \K\ E ( r T 2 ) < C\lt\. (2.25) 

Combining ([2T23]) . (l2?24j) and ([2/25]) . we infer (|272"T|) from (j2?22j) . This, together with 
([2T2Q]) , leads to 



£ r iz,-^_ 1 i 2 dj<c|jt| j 

; = 1 A-l 



which ends the estimate of the second term for the smooth case. 

In general case, let b e , G E , § n ' e , f £ and g e be molifiers of b, G, (j) 71 , / and g, respectively, 
and let (Y e ,Z e ) solution of BDSDE 

Y?* = ^(X^, ...,*£*) +£f(&*)ds + £g e (Zf*)d% - f*Z?*dW s , 0<t<T, 

where X e ' n is the well-know Euler approximation of the diffusion X e , the solution to the cor- 
responding forward SDE (|2.ip modified in an obvious way. Then by the above arguments 
we have 

£ f l \Z* s -Z?^\ 2 ds<C\%\. (2.26) 
Therefore using again Lemma 3.4.2 of Zhang, IfTTTl . we have 

tf \Z S -Zlfds<tf' \Z S -Z?*fds 
< iP[\Z s -ZZf + \ZZ-Z%\ 2 ]ds 

!— 1 Jti-I 



< C<JE [ T \Z s -Z^\ 2 ds+\n\ }. (2.27! 



Applying Lemma 2.2 we have 



limE F \Z s -ZZ\ 2 ds = 0, 



which, combined with (|2.27[) . proves the estimate of the second term of (|3.7|) and together 
with (|2.17p prove the theorem. □ 



3 Discrete-time approximation error 

In order to approximate the solution of the above decoupled FBDSDE (|1.1|) . we introduce 

the following discretized version. Let ji : to < t\ < < t n = T be the partition of the time 

interval [0, T] with mesh 

ItcI = max \tj — 

\<i<n 
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defined in the previous section. Throughout the rest of the paper, we will use the notations. 

Af = ?,-/,■-!, & n Wi = W ti -W ti _ i: and A n Bj = B tj — B tjl for i = 1, ...,n. 
The forward component will be approximated by the classical Euler scheme 



(3.1) 



^0 ~~ ^0! 

K = XL+bplJ AJ+aiXlJAVWi for i=l,...,n 

and we set 

X?=Xl x +b{Xl l ){t-t i ^{)+G{Xl x ){Wt-W ti _ l ) for f e(/,_i,f/). 
We shall denote by {!F™}o<i< n the associated discrete-time filtration define by 

Under the Lipschitz conditions on b and o, the following LP estimate for the error due 
to the Euler scheme is well known 



limsup 1 71 1 ^ 2 max E 



sup \X?-X t \ p + sup \X t -X ti _ u 



p 

KM 'M;_i 

0<?<r f,-i<f<f; 



1/p 

< °°, (3.2) 



for all p > 1 (see e.g Kloeden and Platen, Q). We next consider the following natural 
discrete-time approximation of the backward component Y: 

Y* = h(X?),Z«=0 

2*, =^_ l [(Y t * + g(t u XZ,Y t *)A*B i ) A*Wi], (3.3) 

Y^ = E?_, \Z* + g(?i,X«,Y*) A"*,] +f(t i _ l ,Xl i ,Yl i ,Zl i )Af, (3.4) 

where Ef[.] = E[.|jF f 7t ]. The above conditional expectation are well defined at each step 
of the algorithm. Indeed using the backward induction argument, it easily checked that 
Y t K e L? for all i. 

Remark 3. 1. Using the induction argument, it easily seen that the random variable Y t % and 
Zf are ©i deterministic function of Xf for each i = 0, ...,n. Then using the fixed point of 
Banach argument (|3.4p have a unique solution when the mesh of the partition |tc| is small 
enough. 

For later use, we need a continuous-time approximation of (Y,Z). Since Y t n 
+g(ti,Xj l ,Y t n ) A n Bj = Y t K being in L? for all 1 < i < n, an obvious extension of Ito mar- 
tingale representation theorem yields the existence of the ^-progressively measurable and 
square integrable process Z n satisfying 

f% = E[F* | ] + [' ZfdW s . (3.5) 

Jti-i 
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Then we define inductively 

y* = y^_ { — (t — ti- 1 )f{ti- 1 , Xjf_ [ , ^ j , Zf,_ \)— g{U, X* , ) (B t — B tj , ) 

+ / ZjW,, ^ i < t < U. (3.6) 

The following property of the Z 71 is needed for the proof of the main result of this section. 
Lemma 3.2. For all 1 < i < n, we have 



AfZ^=Et! 



ZJds 



Proof. Since 
recalling (|3.5|> . we have 



A 



-Ef_ , [(y* +g(ti t X* t YZ) A 11 B f ) A" W$] , 



7 71 — _ IF 71 



A*Wi / Zf JW V 



The result follows by Ito's isometry. 

We also need the following estimate, which is a particular case of Lemma 2.4. 
Lemma 3.3. For each 1 <i<n, we define 



□ 



7 K : f 71 



Then 



limsup|7i| ! E 



max sup |y f - Y ti _, \ 2 + Y \Z S -Z% \ 2 ds 

l<'<n;,_,< f <f, i=\ Jt i-i 



< oo. (3.7) 



We are now ready to state our main result of this section, which provides the rate of 
convergence of the approximation scheme (|3.3P and (|3.4|) of the same order than Bouchard 
and Touzi [3]. 



Theorem 3.4. 



Err K (Y,Z) = \ sup E\Y t -Y t K \ z + E 

^0<t<T 



1/2 



< 



ClTll 1 / 2 . 



Proof. In the following, C > will denote the generic constant independent of i and n that 
may take values from line to line. Let i 6 {0, ...,n — 1} be fixed, and set 

8*7, = Y t -Y^,n=Z t -^,^f{t)=f{t,X t J u Z t )-f{t i ,X^,Y^,^) 
and 8"g(0 ^(ajfj-gb^^,), 
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for t £ [tj,tj + i). By Ito's formula, we compute that 

V, = E|5"r f | 2 + Ey^ 1 |5*Z,| 2 dj-|S%. 1 | 2 

= 2E J" + '(5 n Y s ,8 n f(s))ds + J" + '\d K g(s)\ 2 ds, ti<t<t i+l . 

Let P > be a constant to be chosen later. From Lipschitz property of /, g and h, 
together with the inequality ab < fta 2 + b 2 /$ this provides 

v t < ^'' +1 E{[7i| 2 + |x ? -^| 2 + [y 5 -^| 2 + |z 4 .-z^| 2 }^ 



+ J"" CE { |7i| 2 + \X S - X£ +1 1 2 + \Y S - Y£ +l I 2 } ds 



ti+i 



+p / E\8 n Y s \ 2 ds. 



(3.8) 



Now observe that 



E\X s -Xg\ 2 + E\X s -Xg +1 \ 2 <C\%\, 

E\Y S -Y*\ 2 < 2(E|F i -F fi | 2 + E|8 7t F f ,| 2 ) < C (\n\ + E\8%\ 2 ) 
m - Y t f +1 1 2 < 2 (E\Y S - Y ti+l | 2 + E|8 I % +I | 2 ) < C (|jc| +E\5% +1 1 2 ) 



(3.9) 



by (|3 -2 j) and ()3.7j) . Also, with the notation of Lemma 3.3, it follows from Lemma 3.2 that 

E|Z. V -Z*| 2 < 2(E|Z,-Z*| 2 + E|Z*-Z*| 2 ) 



2 E|Z,-Z*| 2 + E 



1 f l M 



E{h n Z r \f ti )dr 



< 



2 (e\Z s -Z*\ 2 + — L ElS^Vj (3.10) 



by Jensen's inequality. 

We now plug (g^J) and (|3~T0|) into (JHH) to obtain 



7 t < E{\k\ + \8%\ 2 + \Z s -Z?f}ds 

+cj t t ' +l E{\K\ + \5%J 2 }ds 

l c r^+i /•'-+• „ , 
+p^ f ' +1 E|8 7t y 4 .| 2 ^ 

< 1 E { |7t| + 16% | 2 + |Z - Z*| 2 } <fo 
+C^'' +1 E{|7t| + |5X 1 | 2 }^ 



C /•'-+' 



+| / E^IVf + p/ E\8*Y s fds. 
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From the definition of V t and (|3.1ip . we see that, for t\ <t < t/ + i, 



E|8 7t F,| 2 + [ t,+1 E\8 % Z s \ 2 ds<$ f t+1 E\8"Y s \ 2 ds+Ai 



(3-11) 



where 



(l+Cn)E\S% +1 \ 2 + ^ 

+- / E\8 n Z s \ 2 ds. 
P Jti 



\n\ 2 + \n\E\Y t n \+ f' +1 E\Z s -Z? 



-ds 



By Gronwall's Lemma, this shows that EjS 71 ^! 2 < A^e^ 71 ' for <t < which plugged 
in the second inequality of ()3.1 1|) provides 



E\S*Y t \' + J E\8 K Z s \ z ds < Ai (I + |jt[p e p W J <A{(1+Cp|jc|) (3.12) 

for 1 71 1 small enough. For t = tj and [3 sufficiently large than C, such that £ < 1, we deduce 
from the last inequality that 

E\8% | 2 + (1 - ~) T +1 E|5*z;| 2 <fe 
P Jti 

< (1+C\%\) U\5% +1 1 2 + |ti| 2 + E[\Z s -Zl\ 2 ]ds\ 

for small |7t|. 

Iterating the last inequality, we get 

E\b*Y ti | 2 + (1 - f' +i E|8*Z,| 2 <fa 
P Ju 

< {l+C\%\) T ^^E\5 n Y T \ 2 + \n\+^^ E[|Z,-Z^J 2 ]rfj|. 



Using the estimate (j3.7|) . together with the Lipschitz property of g and (j3.2|) . this provides 

C rti+i. 



E|8%| 2 + (1-^) f' +1 E\8 n Z s \ 2 ds 
P Ju 

< (1 + C|7c|) I '/ w {E|5*y T | 2 + |n| + C|n| } <C|7t| 
for small |7i|. Summing up inequality f)3- 12|) with t = tj, we get 



(3.13) 



1-^(1+C|3|7T| 

c, 



C E\h n Z s \ 2 ds 
Jo 



< (1 +Cp|jc|) J|jc| + (1 +Cp|7i|)(l +C|7i|)E|5 ,I Fr| 2 



+ 



+ 



(l+cpw^N-i 



E|5%| 2 



(i+cpN)((i+cH) + -|7x|)-i 



n-l 



1=1 



+(1 + CP|7X|)§£ [ ti+1 E\Z s -Z^\ 2 ds. 
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For p sufficiently larger that C, this proves that for small |7i| 



rT 

|2 



Jo 



n-l 

itv |2 



Z,|-V/.s < C ijcl+ElS^rP+W £^(8% 

-l 

+ £E|Z,-Z£| 2 dj 



i=0 

where we recall that C is a generic constant which changes from line to line. We now use 
(f3~T3|) and $SJ$ to see that 

[ T ¥.\§ % Z s \ 2 ds <C\%\. 
Jo 

Together with Lemma [331 and ()3.13|) , this shows that A,- < C|tt| , and therefore, 

sup |5*y ( | 2 <C|ic|. 

o<f<r 

by taking the supremum over t in (13- 12|) . This end the proof of the theorem. □ 
Acknowledgments 

The author would like to thank I. Boufoussi for his valuable comments and suggestions and 
express his deep gratitude to Y. Ouknine and UCAM Mathematics Department for their 
friendly hospitality during my stay in Cadi Ayyad University. 

References 

[1] Bally V., Approximation scheme for solutions of BSDE. Backward stochastic dif- 
ferential equations (Paris, 1995-1996), 177-191, Pitman Res. Notes Math. Ser., 364, 
Longman, Harlow, 1997. 

[2] Bally, V.; Pages, G., A quantization algorithm for solving multi-dimensional discrete- 
time optimal stopping problems. Bernoulli 9 (2003), no. 6, 1003-1049. 

[3] Bouchard B. and Touzi N., Discrete-time approximation and Monte-Carlo simulation 
of backward stochastic differential equations. Stochastic Process. Appl. Ill (2004), 
no. 2, 175-206. 

[4] Douglas, J., Ma, J., Protter, P., Numerical methods for forward-backward stochastic 
differential equations. Ann. Appl. Probab. 6 (1996), no. 3, 940-968. 

[5] Lemor, J.P; Gobet, E.; Warin, X., Rate of convergence of an empirical regression 
method for solving generalized backward stochastic differential equations. Bernoulli 
12 (2006), no. 5, 889-916. 

[6] Karatzas, I. and Shreve, S.E., Brownian motion and stochastic calculus. Graduate 
Texts in Mathematics, 113. Springer-Verlag, New York, 1988. 



16 



[7] Kloeden, RE. and Platen E., Numerical solution of stochastic differential equations. 
Applications of Mathematics (New York), 23. Springer-Verlag, Berlin, 1992. 

[8] Milstein, G.N. and Trekyakov, M.V., Numerical algorithms for forward-backward 
stochastic differential equations. SI AM J. Sci. Comput. 28 (2006), no. 2, 561-582 

[9] Pardoux, E. and Peng, S., Backward doubly stochastic differential equations and sys- 
tems of quasilnear SPDEs Probab. Theory Related Fields. 98 (1994), no. 2, 209-227. 

[10] Yufeng, S.; Weiqiang, Y.; Jing, Y, Numerical computations for backward doubly 
SDEs and , Preprint EXiv:0805.4662 v2 [math.PR] (2008) 

[11] Zhang Jianfeng, Some fine properties of backward stochastic differential equations. 
Ph.D. Thesis, Purdue Univ. (2001). 



17 



